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Recently, spinor Bose-Einstein condensates 
received much attention in both experimental 
BECs have internal degrees of freedom due to the hyperfine spin of t 
a rich variety of phenomena such as spin domains |5| and textures 
valley is so deep that the individual sites are mutually independent, spinor BECs at each 
lattice site behave like spin magnets and can interact with each other through both the 
light-induced and the static, magnetic dipole-dipole interactions. These site-to-site dipolar 
interactions can cause the ferromagnetic phase transition M, M leading to a "macroscopic" 
magnetization of the condensate array and the spin-wave like excitation p, |9| analogous 
to the spin-wave in a ferromagnetic spin chain. For the real spin chain, the site-to-site 
interaction is caused mainly by the exchange interaction, while the dipole-dipole interaction 



is negligibly small. For the spinor BECs in the optical lattice, the exchange interaction 
is absent. The individual spin magnets are coupled by the magnetic and the light-induced 
dipole-dipole interactions js| which are no longer negligible due to the large number of atoms 
N at each lattice site, typically of the order of 1000 or more. Therefore, the spinor BECs 
in an optical lattice offer a totally new environment to study spin dynamics in periodic 
structures. The magnetic soliton excited by the interaction between the spin waves P| is 
an important and interesting phenomenon in spinor BECs. In this paper, we demonstrate 
that the magnetic soliton and elastic soliton collision are admitted for spinor BECs in a 
one-dimensional optical lattice and are controllable by adjusting of the light-induced and 
the magnetic dipole-dipole interactions. 

The Heisenberg model of spin-spin interactions is considered as the starting point for 
understanding many complex magnetic structures in solids. In particular, it explains the 
existence of ferromagnetism and antiferromagnetism at temperatures below the Curie tem- 
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perature. The magnetic soliton 10], which describes localized magnetization, is an important 
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excitation in the Heisenberg spin chain |n|, [TJ, 112, Uj- The Haldane gap of antifer- 
romagnets has been reported in integer Heisenberg spin chain. By means of the neutron 
inelastic scattering [16] and electron spin resonance [3], the magnetic soliton has already 
been probed experimentally in quasi-one dimensional magnetic systems. Solitons can travel 
over long distances with neither attenuation nor change of shape, since the dispersion is com- 
pensated by nonlinear effects. The study of soliton has been conducted in as diverse fields as 
particle physics, molecular biology, geology, oceanography, astrophysics, and nonlinear op- 
tics. Perhaps the most prominent application of solitons is in high-rate telecommunications 
with optical fibers. However, the generation of controllable solitons is an extremely difficult 
task due to the complexity of the conventional magnetic materials. The spinor BECs seems 
an ideal system to serve as a new test ground for studying the nonlinear excitations of spin 
waves both theoretically and experimentally. 

The outline of this paper is organized as follows: In Sec. II the Landau-Lifshitz equation 
of spinor BEC in an optical lattice is derived in detail. Next, we obtain the one-soliton 
solution of spinor BEC in an optical lattice. The result shows that the time-oscillation of the 
amplitude and the size of soliton can be controlled by adjusting of the light-induced dipole- 
dipole interaction. We also present that the magnetization varies with time periodically. In 
Sec. VI, the general two-soliton solution for spinor BEC in an optical lattice is investigated. 



Analysis reveals that elastic soliton collision occurs and there is a phase exchange during 
collision. Finally, our concluding remarks are given in Sec. V. 



II. LANDAU-LIFSHITZ EQUATION OF SPINOR BEC IN AN OPTICAL 



The dynamics of spinor BECs trapped in an optical lattice is primarily governed by three 
types of two-body interactions: spin-dependent collision characterized by the s-wave scat- 
tering length, magnetic dipole-dipole interaction (of the order of Bohr magneton fi B ), and 
light-induced dipole-dipole interaction adjusted by the laser frequency in experiment. Our 
starting point is the Hamiltonian describing an F = 1 spinor condensate at zero tempera- 
ture trapped in an optical lattice, which is subject to the magnetic and the light-induced 
dipole-dipole interactions and is coupled to an external magnetic field via the magnetic 
dipole Hamiltonian H B flQElQ], 



where ip a (r) is the field annihilation operator for an atom in the hyperfine state 
|/ = = a), U B (r) is the lattice potential, the indices a,(3,v,T which run through the 

values —1, 0, 1 denote the Zeeman sublevels of the ground state. The parameter U%£l T (r, r') 
describes the two-body ground-state collisions and U^~p T (r, r') includes the magnetic dipole- 
dipole interaction and the light-induced dipole-dipole interaction. 

When the optical lattice potential is deep enough there is no spatial overlap between 
the condensates at different lattice sites. We can then expand the atomic field operator as 
■0 ( r ) = J2 n E«=o±i ^ol ( n ) n ( r )> where n labels the lattice sites, 4> n (r) is the condensate 
wave function for the nth microtrap and the operators a a {n) satisfy the bosonic commutation 
relations [a a (n) , a/Jl)] = 5 a p&ni- It is assumed that all Zeeman components share the same 
spatial wave function. If the condensates at each lattice site contain the same number of 
atoms N, the ground-state wave functions for different sites have the same form <fr n (r) = 

In this paper we consider a one-dimensional optical lattice along the z-direction, which 



LATTICE 




a 



+ / dvdv'^ a {vW p {v') [ LS*.(r, r') + U^ T (r, r') ] T (r')^(r) + H B , (1) 



a,(3,v,T 



we also choose as the quantization axis. In the absence of spatial overlap between individ- 
ual condensates^and neglecting unimportant constants, we can construct the effective spin 
Hamiltonian M, M as 



H = ~ E A ^ S ™ ' S ' + 2 E X nl^l - ■ B], (2) 

where X n i = 2A^ + A™ z d , A^ d and A^ represent the magnetic and the light-induced dipole- 
dipole interaction respectively. The direction of the magnetic field B is along the one- 
dimensional optical lattice and 7 = gFf^B is the gyromagnetic ratio. The spin operators are 
defined as S n = a]^(n)~F av a v (n) , where F is the vector operator for the hyperfine spin of an 
atom, with components represented by 3 x 3 matrices in the \f = 1, rrif = a) subspace. The 
first term in Eq. (J2J) is resulted from the spin-dependent interatomic collisions at a given 
site, with A^ = (l/2)A a f d 3 r |0 n (r)| 4 , where A a characterizes the spin-dependent s-wave 
collisions. The second and the third terms describe the site to site spin coupling induced by 
the static magnetic field dipolar interaction and the light-induced dipole-dipole interaction. 
For X n i 7^ 0, the transfer of transverse excitation from site to site is allowed, resulting in 
the distortion of the ground-state spin structure. This distortion can propagate and hence 
generate spin waves along the atomic spin chain. For an optical lattice created by blue- 
detuned laser beams, the atoms are trapped in the dark-field nodes of the lattice and the 

n 

light-induced dipole-dipole interaction is very small 8j. However, this small light-induced 
dipole-dipole interaction induces the amplitude and size of the soliton varying with time 
periodically as we will show in the following section. 

From Hamiltonian we can derive the Heisenberg equation of motion at kth site for 
the spin excitations. When the optical lattice is infinitely long and the spin excitations 
are in the long-wavelength limit, i.e., the continuum limit, Sk — > S(z,t), we obtain the 
Landau-Lifshitz equation of a spinor BECs in an optical lattice as follows: 
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where we assume that all nearest-neighbor interactions are the same, namely A <n j> = A, 
which is a good approximation in the one-dimensional optical lattice for the large lattice 



constant Q|. In a rotating frame around z-axis with angular frequency ^ the spin vector 
S is related with the original one by the transformation 

S x = S x ' cos(^t) + S y ' sin(^t), S y = S y> cos(^t) - S x ' sin(^t). (4) 

lb lb IV lb 

Thus the Landau-Lifshitz equation (jSJ) in the rotating coordinate system can be written as 

B B 2 B 2 

—S x = S y —S z - S z —S y - 16p 2 S y S z , 

at Bz 2 Bz 2 

B B 2 B 2 

—S y = S z —S x - S x —S z + 16p 2 S z S x , 

Bt Bz 2 Bz 2 

B B 2 B 2 
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where p 2 = \ ld /(AX), and the superscript ( ' ) is omitted for the sake of pithiness. The 
dimensionless time t and coordinate z in Eq. @ are scaled in unit 2X/h and a respectively, 
where a denotes the lattice constant. Also, the terms including the external magnetic field 
in Eq. (JHJ) have been eliminated with the help of the transformation. 

III. ONE-SOLITON SOLUTION OF SPINOR BEC IN AN OPTICAL LATTICE 

The equation (|SJ) has a form of the Landau-Lifshitz (XL ) ty pe which is similar to the 
LL equation for a spin chain with an easy plane anisotropy 18]. By introducing a partic- 
ular parameter Huang |14| showed that the Jost solutions can be generated and the Lax 
equations are satisfied, and moreover constructed Darboux transformation matrices. An 
explicit expression of the one-soliton solution in terms of elementary functions of z and t 



141 ] we obtain both the one-and 



was reported. Here using the similar method in Ref. |l1 
two-soliton solutions (for detail see the appendix) denoted by S(n) with n = 1, 2 of Eq. (JHJ) 
in the following form: 



Sn = 1 - IT- {X2,n + 2X 3 ,n SU1 2 $„) , 
■'Ml 

S'n = X~ ( X hnVn COSn @n SU1 $ n + X 2 ,n sinn ©n COS $n) , 

S Z n = -T- (Xl,n COsh 9 n COS $ n + X 2 ,n ? 7n sinn ©n SU1 $ n ) , (6) 



where the parameters in the solution are defined by 
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cosh 2 Q n + x 3 n sin 2 $ 
2^4, n [z V n t z n ) , 




(7) 



With «!,„ = /i n (l + p 2 /\C n \ 2 ), «2,n = V n (l ~ p 2 /\( n \ 2 ), K 3 ,n = P„.(l ~ PVlCj 2 ), K 4,n = 



!/„(! + P 2 /|CJ 2 ), Vn = (KJ 2 + P 2 )/(IU 2 - P 2 ), Xl,„ = (2p n ^n) /|C„| 2 , X 2i n = (2^) /|C 



The parameter Vi denotes the velocity of envelope motion of the magnetic soliton. The real 
constants z\ and <f> x represent the center position and the initial phase respectively The 
parameter Ci = Pi + ?zy i is eigenvalue with p l , v\ being the real and imaginary parts. The 
one-soliton solution (JEj) describes a spin precession characterized by four real parameters: 
velocity Vi, phase $1, the center coordinate of the solitary wave z\ and initial phase <p l . 
From the one-soliton solution we obtain the properties of the soliton: (i) both the amplitude 
and the size of the soliton vary with time periodically, as shown in figure 1, in which we 
demonstrate graphically the dynamics of soliton with the parameters chosen as fi l = 0.45 , 
V\ = 0.7, p 2 = 0.375, Z\ = —14 and 4> l = 0.5. This property is resulted from the term p in 
Eq. (JHJ which is determined by the light-induced dipole-dipole interaction. This significant 
observation from the one-soliton solution shows that the time-oscillation of the amplitude 
and the size of soliton can be controlled in practical experiment by adjusting of the light- 
induced dipole-dipole interaction, (ii) the magnetization defined by M 3 = dz(l — Sf) 
varies with time periodically as shown in figure 2. These properties are similar to that 
of the Heisenberg spin chain with an easy plane anisotropy where the dipolar coupling is 
typically several orders of magnitude weaker than the exchange coupling and thus would 
correspond to Curie temperatures much below the observed values. Hence the contribution 
of the dipolar coupling to the spin wave can be neglected in practice. However, for the 
spinor BEC in the optical lattice the exchange interaction is absent and the individual spin 
magnets are coupled by the magnetic and the light-induced dipole-dipole interactions. Due 



and X3n = (,^P 2h, n) /(ICJ 2 — P 2 ) 2 ■ The one-soliton solution, namely S(l) , is simply that 




(8) 



to the large number of atoms iV at each lattice site, these site to site interactions, despite the 
large distance between sites,explain the natural existence of magnetic soliton which agrees 
with the results in Refs. p, u|. 

To see closely the physical significance of one-soliton solution, it is helpful to show the 
parameter-dependence of Euler angles of the magnetization vector which in a spherical 
coordinate is 

Sf(z,t) = cosO, Sf + iSl = wn9exp(i<p). (9) 
From Eqs. © and (JSJ) we find 

|4 + 2x3,1 sin 2 $1 
cosfl = 1 ICl1 



cosh 2 [Fx 1 (z - V x t -zi)]+ x 3 ,i sin 2 $ x ' 

71 

if = — + arctan (r/ 1 tan $ x ) + arctan (tanh 6i) , (10) 

where Fx — 1/(2k 4j1 ) and the maximal amplitude A M = 2 [y\j |d| 2 + I X3 1 1) - When 
|Ci| 2 >> p 2 , the phase (p can be rewritten as 

71 

Lp = — + <j) 1 + k\z — Vt\t + arctan (tanh Gi) , (11) 

where the wave number k x = 2k 31 and the frequency of magnetization precession Qi are 
related by the dispersion law 

fil = h (h + 4p 2 /i 1 / |d I 2 ) - 4« 2)1 K4,i. (12) 
We also see that the position of minimum of energy spectrum fli )TTlin = is located at 

^irnin = ( V /(2 P 2 p 1 ) 2 + 4 K2 , lK4 ,l ICll" " 2?^)/ \Clf ■ 

If amplitude Am approaches zero, namely V\ —>■ 0, the parameter Fx diverges and Eq. 
(HOj) takes the asymptotic form 

71 

cos 9 — > 1, <p — > — + X + fciz — fiit, (13) 

indicating a small linear solution of magnon. In this case the dispersion law reduces to 
Sl x = k 1 +4pVi/ IC1I 2 )- 

IV. ELASTIC SOLITON COLLISION FOR SPINOR BEC IN AN OPTICAL 

LATTICE 



The magnetic soliton collision in spinor BECs is an interesting phenomenon in spin dy- 
namics. Here in terms of a Darboux transformation we first of all give the two-soliton 



solution (for detail see the appendix) of Eq. (J3J) 

S x (2) = S?Sf + R 3 Sf + R 5 S*, 
S y {2) = SfSl + R 4 Sf + ReSf, 

S z (2) = S^S z 2 + R 1 S y 1 +R 2 S z 1 , (14) 

where S*, S% and S* {n = 1,2) are defined in Eq. © and Rj (j = 1,2, ...6) take form as 
follows: 

Ri = -7- (xi 2 cosh 6 2 sinh 9 2 - x 2 2 ?7 2 cos $ 2 sin $ 2 ) , 

R 2 = 1 - ^ (cosh 2 9 2 - sin 2 $ 2 ) , 
A 2 

#3 = -7- (xi 2V2 cosh 6 2 sin $ 2 - x 2 2 sinh 6 2 cos $ 2 ) , 
A 2 

i?4 = 1 - x" [x 2 ,2 sinh 2 2 + (x 2 , 2 + 2x 3 , 2 ) sin 2 $ 2 ] , 

A 2 

^5 = -r-(x 2i2 ?7 2 sin$ 2 sinhe 2 - Xi, 2 cosh6 2 cos$ 2 ), 
A 2 

-Re = -7— (x 2 2 % sin $ 2 cos $ 2 + x x 2 cosh 6 2 sinh 6 2 ), (15) 
A 2 

where 9 2 , $ 2 , fi 2 , A 2 , ?7 2 and x m2 (m = 1,2,3) are defined in Eq. (J7J). The solution (fl4"J) 
describes a general elastic scattering process of two solitary waves with different center 
velocities V\ and V^, different phases <3?i and $ 2 . Before collision, they move towards each 
other, one with velocity V\ and shape variation frequency Q\ and the other with V 2 and 
Q 2 . In order to understand the nature of two-soliton interaction, we analyze the asymptotic 
behavior of two-soliton solution ()14p. Asymptotically, the two-soliton waves ()14|) can be 
written as a combination of two one-soliton waves (JBJ) with different amplitudes and phases. 
The asymptotic form of two-soliton solution in limits t — > — oo and t — > oo is similar to that 
of the one-soliton solution (JBJ). In order to analyze the asymptotic behavior of two-soliton 
solutions (JT4*)) we show first of all the asymptotic behavior of S%, S%, S* (n = 1,2), and 
Rj(j = 1, 2, ...6) in the corresponding limits t — > ±oo from Eqs. (JBJ) and (JTHj) 

-Rl — >■ ±Xl,2J -^2 — * 1 - X2,2> -R3 — >■ 0, 
-R4 — > 1 - X2,2) -^5 — >• 0, -R 6 — >■ TXl,2> 

-»• 1, -> 0, ^ -> 0, as t -> ±00. (16) 

Without loss of generality, we assume that K 4jn > (n = 1, 2) and Vi > V 2 which corresponds 
to a head-on collision of the solitons. For the above parametric choice, the variables Q n (n = 



1,2) for the two-soliton behave asymptotically as (i) Oi ~ 0, G>2 ~ ±00, as t — > ±00; and 
(ii) 2 ~ 0, ©1 ~ Too, as t — > ±00. This leads to the following asymptotic forms for the 
two-soliton solution. (For the other choices of K^ n and V n , similar analysis can be performed 
straightforwardly) . 

(i) Before collision, namely the case of limit t — > — 00. 
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where A = arctan [2/i 2 z/ 2 / (// 2 — ^ 2 )] and the parameters 9 and ip are defined in Eq. (fTTJl) . 



(b) Soliton 2 ( 6 2 ~ 0, a -> 00) 
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ii) After collision, namely the case of limit t - 


-> 00. 



(18) 



(a) Soliton 1 (Q 1 ~ 0, 2 -> 00) 
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(b) Soliton 2 ( 9 2 ~ 0, 0! - 
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here for the expressions of solitons before and after collision, S®, S% and S* (n = 1,2) are 
defined in Eq. (J2J. Analysis reveals that there is no amplitude exchange among three 
components S x , S y and S z for soliton 1 and soliton 2 during collision. However, from Eqs. 
IJI7J) and (fTTJj) one can see that there is a phase exchange 20 A between two components S y 
and S z for soliton 1 during collision. This elastic collision between two magnetic solitons in 
the optical lattice is different from that of coupled nonlinear Schrodinger equations |20j. It 



si 



shows that the information held in each soliton will almost not be disturbed by each other 
in soliton propagation. These properties may have potential application in future quantum 
communication. It should be noted that the inelastic collision may appear if the influence 
of higher-order terms in Eq. (J2J) is considered. 

V. CONCLUSION 

Magnetic soliton dynamics of spinor BECs in an optical lattice is studied in terms of 
a modified Landau-Lifshitz equation which is derived from the effective Hamiltonian of a 
pseudospin chain. The soliton solutions are obtained analytically and the elastic collision 
of two solitons is demonstrated. The significant observation is that time-oscillation of the 
soliton amplitude and size can be controlled by adjusting of the light-induced dipole-dipole 
interactions. 

It should be interesting to discuss how to create the magnetic soliton and how to detect 
such magnetic soliton in experiment. In the previous work p using Landau-Zener rf-sweeps 
at high fields (30 G) a condensate was prepared in the hyperfine state \ f = l,rrif = 0), i.e. 
the the ground state of the spinor BECs. Then the atoms of the ground state can be excited 
to the hyperfine state \f = 1, m/ = ±1) by laser light experimentally. Therefore the excited 
state of the spinor BECs, i.e. the magnetic soliton can be created. As can be seen from 
Fig. 1, the spatial-temporal spin variations in the soliton state are significant. This makes 
it possible to take a direct detection of the magnetic soliton of spinor BECs. By counting 
the difference numbers of the population between the spin +1 and —1 Zeeman sub level, the 
average of spin component < S z > is measured directly. While transverse components can 
be measured by use of a short magnetic pulse to rotate the transverse spin component to 
the longitudinal direction. Any optical or magnetic method which can excite the internal 
transitions between the atomic Zeeman sublevels can be used for this purpose. In current 
experiments in optical lattices, the lattice number is in the range of 10-100, and each lattice 
site can accommodate a few thousand atoms. This leads to a requirement for the frequency 
measurement precision of about 10-100 kHz. This is achievable with current techniques. We 
can also see that the detection of the magnetic soliton of the spinor BECs is different from 
that of the Heisenberg spin chain. 

The magnetic soliton of spinor BECs in an optical lattice is mainly caused by the magnetic 



and the light-induced dipole-dipole interactions between different lattice sites. Since these 
long-range interactions are highly controllable the spinor BECs in optical lattice which is an 
exceedingly clean system can serve as a test ground to study the static and dynamic aspects 
of soliton excitations. 
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VII. APPENDIX 

The corresponding Lax equations for the Eq. © are written as 



d z G (z, t) = LG (z, t) , d t G {z, t) = MG (z, t) , (21) 

where 

L = -ieS z a 3 - iq (S x a l + S y a 2 ) , 

M = i2q 2 S z a 3 + %2qt (S x a 1 + S y a 2 ) - iq {S y d z S z - S z d z S y ) a x 

-iq (S z d z S x - S x d z S z ) a 2 - it {SxdxS" - S y d z S x ) (7 3 , (22) 

here aj(j = 1, 2, 3) is Pauli matrix, the parameters e and q satisfy the relation e 2 = q 2 + 4p 2 . 
Thus Eq. (0) can be recovered from the compatibility condition dtL — d x M + [L, M] = 0. 
We introduce an auxiliary parameter q such that 

e = 2p^— q = 2p r , 23 

q — q q — q 

and the complex parameter is defined by q = (£ + p) / (( — p) . 

It is easily to see that So = (1, 0, 0) is a simplest solution of Eq. (J3J). Under this condition 
the corresponding Jost solution of Eq. (j2~Tj) can be written as 



Go = U exp {—iq (z — 2et) cr 3 } , 



(24) 



where U — \ {I — % io\ + cr 2 + ^3)} with / denoting unit matrix. In the following we con- 
struct the Darboux matrix D n (q) by using the following recursion relation 

G n (q) = D n (q) G n _x (q) , n = 1, 2, 3, (25) 

where D n (g) has poles. Since 



e (~q) = e{q) t q (-q) = -q (q), L (-q) = a x L {q)a 1 , M (-q) = a x M {q)a u (26) 
we then have 

Go (~q) = -ioiGq (q), G n (-q) = -io x G n (q), D n (-q) = a x D n (q)a 1 , (27) 

where the overbar denotes complex conjugate. Suppose that q n is a simple pole of D Tl (q) , 
then — q n is also a pole of D n (q). If D n (q) has only these two simple poles we have 

D n (q) = C n P n (q) , (28) 

Pn (q) = 1 + ^^Pn + q -^-Pn, (29) 

q - q n q + q n 

where C n , P n , and P n are 2x2 matrix independent of q, the terms (q n — q n ) C n P n and 
(q n — q n ) C n P n are residues at q n and q n , respectively. From Eq (}2*Tj) . we have 

C n = a 1 C n a 1 ,P n = (J\P n (J\. (30) 

From Eqs. (J22J) and (23) we see that 

L(q) = -Lt (q) , M (q) = — (q) , G \q) = G\ (q) , (31) 

and hence we have 

G-\q) = G{ (q) , D-\q) = £>t (3) = pt (g) tft. ( 32 ) 
Since .D n (g) I?" 1 (q) = D n (q) D^ 1 (q) = I, it has no poles. Then we obtain 

PnPl (q n ) = 0, P n (i - Pi + ^q^Pl) = 0, (33) 



which shows the degeneracy of P n . One can write P n = ( g n w n ) T ( T ra £ n ) where the 
superscript T means transpose. Substituting this expression into (jHHJ) we obtain 



Pn(q) 



q n \ T n\ 2 + q n \U 2 o 



A n (q-q n )(q + q n ) ^ ?n |T„| 2 + ? n |£ n | 2 

x I 2 | ^ n l^nl J _|_ g ^g 2 —2^ I ^ f «£r 

g„ |T„| 2 + g„ |^„,| 2 / \ £ n T n 

|?n| 2 ( gJej2+ ^ |T " 12 ^ | (34) 

gn|T„| 2 + gj£ n | 2 



where 



An = kn| 2 (|T n | 2 + |^ n | 2 ) 2 + \q n - q n \ 2 |T n | 2 |£ n | 2 . (35) 

To determine £ n and T n , we substitute (J25J) into (f2"TJ) and take the limit q — ► g n and then 
obtain 

<9 z L> n (g) = L n (q) D n (q) - D n (q) L n _i (q) , 

«9 t D n (q) = M n (q) D n (q) - D n (q) M n ^ (q) . (36) 

Because of the degeneracy of P n , the second factor of the right-hand sides of Eq. (jSEJ), 
namely, ( T n £ n )G n -i (q n ) must appear in the left-hand side in its original form and, hence, 
it is independent of z and t. We simply let 

( T n i n ) = ( b n 1 ) G-\ (q n ) , (37) 

here b n is a constant. Hence, the Darboux matrices D n (q) have been determined recursively, 
except for C n . By a simple algebraic procedure, it is seen that A n is always non-vanishing 
regardless of the values z and t. This shows the regularity of P n and then P n (q). In the 
limit as q — > 1, from Eq. (|23~j) we have 

e(q), q (q)^2p^— + 0(l), 
q - 1 

and then from Eq. (J3fij) we obtain 

S (n) ■ a = D n (1) [3 (n - 1) • a] D\ (1) , n = 1, 2, 3, .... (38) 
Considering the Eq. ()31)1 and Eq. ()32)1 we get 

C n Ct = I, (39) 



which shows that the matrix C n is diagonal with the help of the Eq. (JHUj) and 

(C n ) u =(C- n ) 22 ,\(C n ) n \ = l, (40) 

then we can write C n = exp (iuj n a^/2) which is real and characterizes the rotation-angle of 
spin in the iq/-plane. It is necessary to mention that uo n may he dependent on z and t. To 
determine u n one must examine the Lax equations carefully. Since exp (iu n 03/2) denotes 
a rotation around the z-axis, it does not affect the value of S z . Substituting ()28j) into (JHHJ) 
and taking the limits as q — > 00 and q — > respectively, we obtain 

d z {C n P n (0)} = L n (q) {C n P n (0)} - {C n P n (0)} L n _! (q) , 

d z {C n } = -i2pS z (n)a 3 {C n } + {C n }i2pS z (n-l)a 3 . 

Comparing these two equations, we derive 

C^l^'l^ ° ). (4!) 

The Eq. (jSU) gives 

T n = f n + i/ n , £ n = / n — i/„ , (42) 

where 

= exp (-6„ + l$n) , 

here the parameters n and $ n are defined in Eq. 0. Setting n = 1 and substituting 
the Eqs. flUD, JSU), (JUJ), and into (JSHI) we can obtain the one-soliton solution ©. 
Setting n = 2 and with the similar procedure the expression of the two-soliton solution (JUjl 
is derived. 
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Figure Captions 



Figure 1 The amplitude and size of soliton in Eq. (jHJ) vary periodically with time, where 
p 1 = 0.45, V\ = 0.7, z\ = —14, (p 1 = 0.5, p 2 = 0.375. The unit for time t is 2X/H and a for 
space z. 

Figure 2 The magnetization M3 (the integral dz(l—Sf)) vary with time periodically, 
where fx 1 = 0.45, V\ = 0.7, z\ = —14, 1 = 0.5, p 2 = 0.375. The unit for time t is 2X/H and 
a for space z. 
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